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Abstract 

Short-wave perturbations in a relaxing medium, governed by a special reduction of the Ostro- 
vsky evolution equation, and later derived by Whitham, are studied using the gradient-holonomic 
integrability algorithm. The bi-Hamiltonicity and complete integrability of the corresponding dy- 
namical system is stated and an infinite hierarchy of commuting to each other conservation laws 
of dispersive type are found. The two- and four-dimensional invariant reductions are studied in 
detail. The well defined regularization of the model is constructed and its Lax type integrability 
is discussed. 
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1. Introduction 

Many important problems of propagating waves in nonlinear media with distributed parameters 
can be described by means of evolution differential equations of special type. In particular, if 
the nonlinear medium is endowed still with spatial memory properties, the propagation of the 
corresponding waves can be modeled by means of the so called generalized Ostrovsky evolution 
equations [1]. It is also well known [5, 7, 6] that shortwave perturbations in a relaxing one 
dimensional medium can be described by means of some reduction of the Ostrovsky equations, 
coinciding with the Whitham type evolution equation 

(1.1) du j ' dt = 2uu x + / K.(x,s)u s ds, 

Jr 

discussed first in [5]. Here the kernel K, : R x R — > R depends on the medium elasticity properties 
with spatial memory and can, in general, be a function of the pressure gradient u x G C°°(R;R), 
evolving in respect to equation (1.1). In particular, if K.(x, s) = \ \ x — s |, x, s € R, then equation 

(1.1) can be reduced to 

(1.2) du/dt = 2uu x + d~ V 

which was, in particular, studied before in [3, 2, 6, 7, 5]. 

Since some media possess elasticity properties depending strongly on the spatial pressure gra- 
dient u x , x G R, the corresponding Whitham kernel looks like 

(1.3) JC(x, s) := — 9(x — s)u s 

for i,s£R, naturally modeling the relaxing spatial memory effects. The resulting equation (1.1) 
with the kernel (1.3) becomes 

(1.4) du/dt = 2uu x - d~ x u\ := K[u], 

which appears to possess very interesting mathematical properties. The latter will be the main 
topic of the next sections below. 

2. Lagrangian analysis 
A more mathematically correct form of equation (1.4) looks like 

(2.1) u xt = 2(uu x ) x - ul, 

being a nonlinear hyperbolic flow on the axis R. Concerning the preceding form of (1.4) it is 
necessary to define the operation d^ 1 : C°°(]K;M) — > C°°(R;R), which is not an easy problem. 
Since equation (2.1) is well defined in the space of 2n- periodic functions C^RjR), one can 
determine upon its subspace C^(R;R) C Co^(R;R) of functions under the condition f Q f(s)ds 
= for any / G (^(RjR) the inverse operation 

(2.2) gr\.):=-[^ (.)ds - £ (.)*], 
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being definite for all x € R and satisfying the natural property d ■ d^ 1 = 1. Thereby, for 
convenience, we will consider the flow (1.4) as that in the smooth functional submanifold 
M := C^(R;R). The corresponding vector field K : M -> T(M) defines on M a dynamical 
system, which appears to possess both Lagrangian and Hamiltonian properties. 

To demonstrate them in detail, consider the partial differential equation (2.1) and prove that 
it is of Lagrangian form, that is 

m 

(2.3) Uxt = —^:=Z[u], 

where H$ : M — > R is some Frechet smooth Lagrangian function. To prove (2.3), following the 
scheme in [8, 14, 12], it is enough to state that only the Volterrian identity £' = £'* holds, that is 

(2.4) [2(uu x ) x - u 2 J = [2(uu x ) x - u 2 J*, 

where the sign ""' means the Frechet derivative with respect to the variable u G M and " * " 
means the corresponding conjugation with respect to the natural scalar product on the tangent 
space T(M) ~ T*(M). As a result, there exists a Lagrangian function H$ : M — ► R in the 
following explicit form: 

(•lit r27T 

(2.5) Htf := / Tifidx = / uu x dx, 

Jo Jo 

where we used the standard [8, 14] homotopy formula H$ = Q d\{gradH$[uX) ,u). Thus, expres- 
sion (2.3) can be presented as the Euler equation 

(2.6) 6£/6u = 0, 
where, by definition, 

(2.7) C := / / {-u x u T - Hn)dxdT. 

Jo Jo 2 

Recall now, that owing to the standard results [8, 12, 14, 19, 20], any Lagrangian system in the 
form (2.6) is Hamiltonian. To show this, rewrite the action functional (2.7) as 

(2.8) C= f [(<p,u T ) - H#}dr, 

Jo 

where (p := (1/2)^ G T*(M). Then the condition (2.6) gives rise to the equality 

(2.9) ut = grad H#[u] = K[u], 
where, by definition, 

(2.10) iT 1 := (ff - if'* = d/dx. 

As it is easy to see, the operator •& := d^ 1 : T*(M) — ► T(M) is necessary implectic [8, 10, 12] and 
with respect to the flow (2.9) also Noetherian. Thus we have stated the following [16] theorem. 

Theorem 2.1. The partial differential equation (2.1) is equivalent on the functional manifold 
M to the Hamiltonian flow (2.9) with the Hamiltonian function (2.5) and co-implectic structure 
(2.10). 



This result means that our flow (1.4) on M, being Hamiltonian, is conservative, thereby one 
can expect it possesses also an additional hidden infinite hierarchy of conservation laws, which is 
very important [14, 8, 10, 18, 11] for its integrability analysis. This assumption, as we shall show 
below, appears to hold really. 

3. Gradient-holonomic analysis 
Since any conservation law 7 G D(M) satisfies the linear Lax equation 

(3.1) dil>/dt + K'*il> = 0, 

where ip = grad 7 € T* (M) , under the condition of its existence in the form of a local functional 
on M, it can be found for instance, by means of the asymptotic small parameter method [8]. In 
particular, one easily gets that expressions 

(3.2) fo = U xx , Ipr,-! = g(«x - («)L) 

satisfy the Lax equation (3.1) and are the gradients of the corresponding functionals on M, that 
is 

(3.3) ip# = grad 70 Wi = grad 7^, 
where 

(3.4) 7,? = - J u\dx 7^ = 2 _/ o uu l dx - 

Thus, we have stated that our dynamical system (1.4) allows additional invariants (conservation 
laws), which can be used within the gradient-holonomic algorithm [8, 15, 18] for finding new as- 
sociated nontrivial implectic structures on the manifold M. Namely, let us represent conservation 
laws (3.2) in the scalar product form on M as 

(3.5) 70 = (W,«x) 7t7-i = (tfri-i^x), 

where 

(3.6) <P# = \u x , <p r ,_ l =-±dr 1 u 2 x €T*(M). 
Then operators 

v- 1 = = \d-{-\d) = d, 

(3.7) r]Z\ = - ^*_! = d~ X u xx + u xx d~ x 

will be co-implectic [8, 12, 14] on M, and, as it is easy to check, also Noetherian with respect to our 

dynamical system (1.4). Moreover, via direct calculations one can show that the corresponding 

implectic operators 1), rj-i : T*(M) — > T{M) are compatible on M, that is for any A £ R the 

expression $ + Xr/-i is implectic too on M [8, 9, 12]. Really, it is enough to show [12, 8] that 

the operator "!? _1 r/_i^ _1 : T*(M) — > T(M) is symplectic on M, that is the differential two-form 
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Q( 2 ) : = J 27r dx(du A d- 1 ri_ l 'd- l du) G A 2 (M) is closed, or dQ^ = 0. The latter equality is easily 
checked by direct calculations. This means, in particular, that all operators of the form 

(3.8) Vn = $( V ZI$T 

for n G Z will also be implectic on M. Another consequence from this fact is the existence of an 
infinite hierarchy of invariants 7„ G D(M), neZ, satisfying the expressions 

(3.9) K[u] = -r] n grad -y n . 

As a particular case one can define an implectic operator n : T*(M) — ► T(M) in the form 

(3.10) r? = tfrCitf = d-^d- 1 + cT 1 ^- 2 . 
Whence from (2.9) we obtain that 

(3.11) ut = K[u] = —figrad H$ = —rjgrad H v , 
where 

H# = / uu 2 x dx, H v = / u 2 x dx. 
Jo Jo 

The set of expressions (3.8) can be equivalently rewritten in another useful form as 

(3.12) Ai? grad 7(A) = 77 grad 7(A), 

being in some sense equivalent [8, 14, 12] together with equation (3.1) to the adjoint Lax type 
representation 

(3.13) dA/dt = [A, K'*] 

for the dynamical system (1.4), where A := i? _1 r/ : T*(M) -> T*(M) is a so called [14, 12, 8, 12] 
recursion operator and 7(A) G D(M), A G C, is a generating function of the infinite hierarchy of 
conservation laws (1.4). In particular, as |A| — > 00 the asymptotic expansion 

(3.14) grad 7(A)| |AHoo - ^ 5rad 7 i 
holds, where 

(3.15) grad 7„ = A n ^rad 70, 70 := H v , 

for all n G Z+. Concerning this infinite hierarchy of conservation laws one can easily check, that 
all of them are dispersiveless. The result obtained above can be formulated as the next theorem. 

Theorem 3.1. The dynamical system (1-4) on the functional manifold M is a compatible bi- 
Hamiltonian flow, possessing an infinite hierarchy of commuting functionally independent dis- 
persionless conservation laws, satisfying the fundamental gradient identity (3.12). The latter is 
equivalent together with the relationship (3.1) to the adjoint Lax type representation (3.13). 

As was mentioned above, the hierarchy of commuting flows K n := —1? grad j n , n G Z + , shows 

an interesting property of their dispersionless. In particular, this entails that they can not be 
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treated effectively by means of the gradient-holonomic algorithm [8, 15, 18]. In particular, the 
corresponding asymptotic solutions to the Lax equations 

(3.16) dcp/dT n + = 0, iff + <p' + , 

where | A |— > oo and du/dr n = K n [u], t„ £ I, n £ Z, do not give rise to explicit functional 
expressions, defining a new associated hierarchy of conservation laws for the dynamical system 
(1.4). Nonetheless, the corresponding hierarchy of dispersive commuting flows on M does exist for 
(1.4), being simply associated with the trivial flow du/dto := on M. Namely, let Hq G D(M) be 
a conservation law of (1.4), satisfying the kernel condition for the operator rj : T*(M) — ► T(M), 
that is 

(3.17) du/dto = 0:=rj grad H . 

It is easy to find from (3.17) and (3.10) that grad H = ^{u^)- 1 / 2 ]^ G T*(M), whence 

/•27T 

(3.18) H = 4 Ju— X dx. 

Jo 

The obtained invariant (3.18) allows to construct a new associated with (1.4) commuting flow 

(3.19) du/dr = -& grad H = u xxx (u xx )~ 3/2 := K[u], 

t G E, which, as it is easy to see, already possesses a nontrivial dispersion. This means that the 
Lax equation 

(3.20) dip/dr + k'*ip = 0, 

allows as |A| —>■ co an asymptotic solution 92 := ip(r, x; A) G T*(M) (g) C, where 

(3.21) <p(t,x;\) ~ exp(X 3 T+ I a(y;X)dy), 

Jxo 

a(x-X) ~ ^2 <7j-i[u]\~ j+1 . 

The nontrivial functionals 7j_i := J 2w aj-i[u]dx, j G Z + , are, obviously, functionally indepen- 
dent and commuting conservation laws both of the dynamical system (3.19) and of our dynamical 
system (1.4). As a result of some simple but slightly tedious calculations one finds that 

(3-22) cr_i = cr = ^u~lu xxx , a x = ^{u xx )~ b/2 u 2 xxx , 

and the corresponding hierarchy of already dispersive invariants is given as 

/•27T 

(3.23) 7_i = / dx > 7o = 0, 

J 

/•27T 

71 = I / «^ /2 «xxa; 72 = 0, ... , 

JO 

and so on. Then, owing to conditions (3.19) and (3.20) the generating functional 7(A) := 
J Q 2w a(x; X)dx, X G C, satisfies [8, 12, 15] the following gradient relationship 

(3.24) A 2 i? grad 7(A) = 77 c/md 7(A), 
suitably modifying the relationship (3.12). 



The obtained results are very important for further analytical studying Lax type integrability 
of the dynamical system (1.4) and finding, in particular, a wide class of its special soliton like 
and quasi-periodic solutions by means of analytical quadratures. Some of these aspects of the 
integrability problem are presented in the section below. 

4. Lax type representation and finite dimensional reductions 

Since the functional solution (3.21) satisfies the Lax type equation (3.20), it can be considered 
[15, 18, 10] as a Bloch type eigenfunction of the adjoint Lax type representation (3.13), that is 

(4.1) A<p(x;X) = \ 2 <p(x;\) 

for all A G C and i£R. This gives rise, following the gradient-holonomic algorithm [15, 18], to the 

existence of a standard Lax type representation for the associated dynamical system (3.19) and, 
thereby, for our Whitham type dynamical system (1.4). Omitting here the related calculations, 
we find surprisingly that this adjoint to (4.1) Lax type spectral problem for the flow (1.4) is 
equal to 

(4.2) Lf := ( V ) /, 

where an eigenfunction / 6 Loo(IR;C 2 ) and A € C is a time independent spectral parameter. The 
result (4.2) can be effectively used for solving our nonlinear equation (1.4), making use either of 
the inverse spectral transform method [18, 17, 15, 11, 10] or of the dual Bogoyavlensky-Novikov 
method [18, 8] of finite dimensional reductions. For the latter case we need to construct an finite 
dimensional invariant symplectic functional submanifolds M 2N C M, N € Z+, and to represent 
the main vector fields d/dx and d/dt on them as the corresponding commuting to each other 

Hamiltonian flows. Moreover, since these flows on M 2N appear to be Liouville- Arnold integrable, 
we obtain both the complete integrability of our dynamical system (1.4) in quadratures and their 
exact solutions, expressed, in general, by means of Riemannian theta-functions [18, 17, 15, 10] on 
some specially constructed algebraic Riemannian surfaces. 

Below we consider, for simplicity, the following invariant two- and four-dimensional functional 
submanifolds: 

(4.3) i) M 2 := [tie M : gradC 2 [u] = 0}, 
where £2 : = H$ + c v H n £ V(M), and 

(4.4) it) M 4 := {u G M : gradC A [u] = 0}, 

where £4 := 7_i + c#H$ + CqH^ 6 V(M). 

Case i). We have, therefore [8, 13, 18, 17], on the invariant manifold M 2 commuting Hamil- 
tonian vector fields d/dx and d/dt with respect to the canonical symplectic structure 



(4.5) 



( 2 ) :=daW, 
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where 1-form G A 1 (M) is determined by the Gelfand-Dickey [13, 8] relationship 

(4.6) CLC2M = gradC2[u]du + da^ /dx, 
holding on M. One now easily finds, that for all u G M 2 C M 

(4.7) grad£ 2 [u] = u 2 , - 2(u x u) x - 2c v u xx = 0, 

= 2(u + c v )u x du, 
a/ 2 -* = d[2(u + c v )u x ] A (in := dp A c/g, 

where we have put p := 2(u + c^)u x and g := u. The corresponding Hamiltonian functions 
and M*) G V(M 2 ) for Hamiltonian flows 

(4.8) dq/dx = dh^/dp, dp/dx = -dh^ /dq, 

dq/dt = dh^/dp, dp/dt = -dh^/dq 

are found [8], respectively, from the determining relationships 

(4.9) grad£ 2 [u] u x := -dh^ /dx, gradC 2 [u] u t := -dh^ /dx, 
whence we get that 

(4.10) h& = (u + c v )u x = A( P [ h® = -2c v (u + c v )u x = 7 C f 2 . 

4{q + c v ) 2(q + c v ) 

One sees easily that two flows d/dx and d/dt on the two-dimensional invariant submanifold M 2 C 
M of infinite period are proportional, confirming the classical fact [20, 19]: upon two-dimensional 
symplectic manifold there exist the only functionally-independent invariant commuting to each 
other. 

The set of Hamiltonian equations (4.9) for the flows d/dx and d/dt has the simple form 
(4.11) 



dq 


p 


dp 


p 2 


dx 


2(9 + Cr,) ' 


dx 


<q + c v ) 2 


dq 


-c v p 


dp 


-c v p 2 


dt 


(q + c v ) ' 


dt ~ 


2(q + c v ) 2 > 



whose solution, given by the exact formula 
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(4.12) q(x,t) =-c + [-VhW(x - 2c v t) + k} 2/3 u(x,t), 

with k £ H being some real constant, supplies us, evidently, with an exact partial one-parametric 
solution to our Whitham type nonlinear equation (1.4). 

Case ii). Similarly as above, we find the quantities 

(4.13) gradC 4 [u] = ( ) xx + c#[ul - 2(uu x ) x ] - 2c v u xx = 0, 

= [2(c#u + c v )u x - {—^—) x ]du + 1 — du x , 



whence the symplectic structure is given as 

(4.14) J 2) = d[2(c^u + c v )u x - {——=) x \ Adu + d( - ) A du x 

= dpi A dqi + dp2 A dq 2 , 

where we put, by definition, (qi := u, q 2 : u x , p\ := 2 (c#u + c^)^ - ( 2 J Uxx ) x , P2 ■= 2 J Uxx ) £ 
which are canonical symplectic coordinates on the invariant functional submanifold M 4 C M. The 
commuting to each other Hamiltonian functions, related with flows d/dx and d/dt, are equal to 
the next algebraic expressions 

(4.15) = ul(c#u + c v ) - ( - — ) x u x - y/u^ 

= ^Q 2 {cMi + c v ) - q 2 pi - 1/ {2p 2 ) 

and 

(4.16) = qi/(2p 2 ) - 2q 1 q 2 \p 1 - q 2 {c$qi + c v )}. 

As a result, we have reduced our Whitham type dynamical system (1.4) upon the constructed 
four-dimensional invariant submanifold M 4 C M, on which it is exactly equivalent to two com- 
muting canonical Hamiltonian flows 

dqj/dx = dhW/dpj, dpj/dx = -dhW/dqj 7 

(4.17) dqj/dt = dh^/dpj, dpj/dt = -dh^/dqj 

for j = 1, 2, where the corresponding Poisson bracket {h^ x \ h®} = on M 4 . Thereby, owing to 
the classical Liouville- Arnold theorem [20, 19, 8] our Whitham type dynamical system (1.4), 
reduced invariantly upon the four-dimensional invariant submanifold M 4 C M, is completely 
integrable by quadratures. This result we will formulate as a final theorem. 

Theorem 4.1. The Whitham type dynamical system (1-4), reduced upon the invariant two- 
parametric four- dimensional functional submanifold M 4 C M is exactly equivalent to the set of 
two commuting to each other canonical Hamiltonian flows (4-17), which are completely Liouville- 
Arnold integrable by quadratures systems. The corresponding Hamiltonian functions are given by 
expressions (4-15) and (4-16). 

The results obtained above make it possible to construct a wide class of exact two-parametric 
solutions of the Whitham type nonlinear equation (1.4) by means of quadratures. This very 
interesting and important problem we plan to investigate in detail in another paper. 

5. REGULARIZATION SCHEME AND THE INTEGRABILITY PROBLEM 

Define a smooth periodic unction v G C^(R;M), such that 
(5.1) v := d- x u\ 

for any i,tel, where the function u 6 C|£ (R; R) solves equation (1.1). Then it is easy to state 

9 



that the following nonlinear dynamical system 

/r n\ Ut = 2UU X — V ] r ,r ! 

of hydrodynamic type, being already well defined on the extended 27r-periodic functional space 
M := Cg^EjE), is completely equivalent to that given by expression (1.1). Thereby, mapping 
(5.1) regularizes the previously not completely determined expression (1.1), making it possible to 
pose a new integrability problem for dispersionless dynamical system (5.2) of hydrodynamic type 
in the functional space M. To proceed, one can analyze this problem by means of the gradient- 
holonomic method [8, 15], which is similar to the approach applied above for studying dynamical 
system (1.1). Based on preliminary performed calculations and obtained analytical properties of 
dynamical systems (5.2) we can state that it is also an integrable flow in the functional space 
M, possessing a suitable Lax type representation. We plan to investigate this problem in detail 
in another work under preparation. 
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